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■ A short presentation of the results of the analytical evaluation of the O(a^) QCD 

. contributions in the fermionic decay rates of the neutral CP-odd minimal super symmetric 

I Higgs particle is made. The corrections due to the nonvanishing quark masses are 

' included. The results are presented both in terms of running and pole quark masses. 
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A crucial test of the most attractive extension of the Standard Model (SM) - the Minimal 
Supersymmetric Standard Model (MSSM) can possibly be done at LEP200. Indeed, as was 
shown in there is a significant mass parameter space increase in the two Higgs doublet 
I model at LEP200, while for the SM Higgs the mass range will increase only marginally. For 
the MSSM, a large theoretically allowed parameter region can be covered 0] with ^/s = 210 
GeV, highest expected luminosity of 500 pb~^ and with the latest CDF result for the top mass 
around 174 GeV (For earlier analyses of the MSSM Higgs discovery potential at LEP and 
LHC see and references therein.) 

To study the detectability of various Higgs particles, it is important to start with the 
search of qualitative as well as quantitive information about their decay processes. One of the 
most interesting experimentally accessible quantities - the branching fractions of the quark- 
antiquark decay modes can be calculated within the perturbation theory. (For a review of the 
SM Higgs phenomenology, taking into account radiative effects see and references therein.) 
A precise theoretical evaluation of decay rates is particularly important from a viewpoint of 
distinction of MSSM and minimal SM Higgs signals. In the recent work [Q the fermionic decay 
rates of the minimal SM Higgs particle have been evaluated to the O(a^) in perturbative QCD. 
In the present work, using the same technique, the 0{a1) QCD contributions to the decay 



rates of the MSSM pseudoscalar Higgs particle into the quark-antiquark pairs are evaluated 
analytically, taking into account the nonvanishing quark masses. The assumption was made, 
that all superpartners are too heavy to appear in the decays of the Higgs boson. 

The minimal supersymmetric extension of the SM is the two-Higgs doublet model. The 
physical Higgs boson spectrum of the MSSM consists of five states: a charged Higgs pair H"^, 
neutral CP-even scalars Hg, H'^ and a neutral CP-odd pseudoscalar Hp (for a review see, e.g., 
the textbook |^). The Lagrangian density, describing the Yukawa interaction of the neutral 
Higgs boson with quarks and leptons has the following form: 

L = {V2GpY'^ T.{CH.ff{a, /3)mfqfqjH, + CH,fj{P)mfqjtj,qjHp) 
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Below qf denote the quark fields with flavor / and mass nif. The coefficients C//^//(a, /3) 
and CHpffiP) are the functions of the ratio of the Higgs- field vacuum expectation values - 
tan/5 and a mixing angle a in the neutral CP-even sector. The expressions for Cn^ffici, (3) 
and CHpffiP) for the particular quark or lepton pair can be found, e.g., in 0. In the SM 
limit Cif^//(a, P) = —1 and the pseudoscalar term drops out in the eq.(0). The corresponding 
calculation of the 0{a'^) SM Higgs decay rate (relevant for the MSSM Higgs scalar as well) 
has been done in p. In the present work one considers the pseudoscalar Higgs, expecting 
that the reader is familiar with the ref. 0. 

The two-point correlation function of the pseudoscalar currents jj^'' has the following form: 
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For the decay width one has: 
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The full 0{as) analytical result for the decay rate of Hp qjqj- in terms of pole quark 
masses looks like 1711: 
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where: 
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The expansion of the r.h.s of eq.(|^) in a power series in terms of small m'j/Mjj^ has the 
following form: 
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where the period covers high order terms ~ {rrif / M^^Y^ , k = 2, 3.... 

The evaluation of the O(a^) corrections to the eqs. (^,||) has been done in a close analogy 
of the similar calculation of the decay rate of SM Higgs boson |^. Namely, the full two- 
point correlation function was expanded in powers of m'j/Q'^ in the "deep" Euclidean 
region (Q^ ^ m'j). The first two coefficient functions in this expansion was evaluated to 
0{al) by applying an appropriate projectors to the Il{Q'^,mf ,m^). is the mass of the 
quark appearing virtually in some topological types of three- loop diagrams and "B" labels the 
unrenormalized quantities. In the calculation of divergent Feynman integrals the dimensional 
regularization formalism |^ and the minimal subtraction prescription |^ in its modified form 
are used. The 75 matrix is defined in D-dimensional space-time as an object with 

{75,7m} = 0' 7575 = 1. 

The above definition causes no problems in dimensional regularization when there are two 75 
matrices in a closed fermionic loop. The corresponding one-, two- and three-loop diagrams 
have been evaluated analytically, using the special computer program HEPLoops The 
graph-by-graph results were summed up with an appropriate symmetry and gauge group 
weights. In the obtained expressions for the Ilj in terms of bare quantities one renormalizes the 
coupling and the quark mass in a standard way (see, e.g., p). Finally, one should analytically 
continue the obtained Ilj from Euclidean to Minkowski space and take the imaginary part at 
s = Mfj^ ( eq.@ ). One obtains the following analytical result for the standard QCD with 
SUc{S) gauge group: 
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where the Riemann functions ({2) = and ({3) = 1.202056903. The last term in eq.(^ 
represents the contributions from the three-loop diagrams containing the virtual quark loop 
(see fig.l of the ref.^). The "triangle anomaly" type correction (see fig.2 of the ref.^) 
requires a special treatment because of 75 in dimensional regularization, and is not included 
here. However, it is reasonable to expect that their contribution will not exceed 1%. Note 
that the "triangle anomaly" type contributions vanish identically in the zero quark mass limit. 
Note also the smallness of corresponding corrections in the scalar channel 

The leading three-loop term ~ mi (the massless approximation) coincides with the one 



obtained in 12, 13 



while the three- loop results ~ m'f,m'jml are new. 



The calculated decay width obeys the homogeneous renormalization group equation: 
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For the standard definition of the QCD /3-function and the mass anomalous dimension see the, 
e.g., ref.|^ (eqs. (14)-(17)). The renormalization group relates the coefficients of log terms in 
eq.(^). The corresponding relations up to O(a^) for the SM Higgs decay, which are valid for 
the pseudoscalar Higgs as well, has been obtained in . 

The solution of the renormalization group equation (0) at fij^ = Mfj^ has the following 
form: 
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For the standard parametrization of the running couphng and the running mass see, e.g., 
ref.p]. One can see that the three-loop mass correction is positive in contrary to the analogous 
result for the SM Higgs p. 

The relation between the MS quark mass mf{M) renormalized at arbitrary M and evalu- 
ated for the A^- flavor theory and the pole quark mass nif has the following form |^: 
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where (3q^\ 7o and are the perturbative coefficients of the QCD (3 function and the mass 
anomalous dimension 7^ correspondingly and can be found, e.g., in ref. Kj and 6{mf,mf') 

see [|): 



can be obtained from the on-shell results of 1114 
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L±{r) can be evaluated for different quark mass ratios r numerically (see the table 1 of the 
ref.0). In the above equations the number of participating quark flavors is specified 
according to the size of M and has no correlation with the quark mass ruf. 
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Substituting eqs. at M = Mh^ and the coefficients of the renormahzation group 

functions into the eq.(P), one obtains the general form for the decay rate THp^qf-q^ in terms of 
the pole quark masses: 
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while the 



The above result con&ms the asymptotic form (|^) of the two-loop exact result ( 
0(0^) expression is new. 

The numerical values of A(m;/mj), defined in the eq.(0), are given in the table 1 of the 
ref.0. The quark masses are estimated from QCD sum rules: rrif, = 4.72 GeV [1^, rric = 1.46 
GeV |T§, = 0.27 GeV |71,_m„ + « m,/13 = 0.02 GeV [0. 

For the decay mode Hp bb with five participating quark flavors one obtains numerically: 
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where 5'-^^ is defined in eq.(^). 

In ref . it was observed that the high order corrections reduce the scale dependence signif- 
icantly and resolve the large discrepancy between the results for the decay rate of SM Higgs 
in terms of running and pole quark masses. The theoretical uncertainty was estimated at 
5%. The same effects were observed for the decay rate of pseudoscalar Higgs particle. In the 
present calculation one also estimates the theoretical uncertainty at approximately 5%. 

It should be stressed, that the virtual top quark, which appears in some topological types 
of three-loop diagrams (see fig.l and fig.2 in ref.|P) may give a nonnegligible contribution. In 
addition, one may include a heavy virtual superpartners. However, the superpartners decouple 
I9| if the Higgs mass is much smaller than the supersymmetry scale. 



Acknowledgments It is a pleasure to thank D.Soper for discussions, D.Broadhurst and 
B.Kniehl for helpful communications. This work was supported by the U.S. Department of 
Energy under grant No. DE-FG06-85ER-40224. 



References 

[1] A. Sopczak, Intern. J. Mod. Phys. A9 (1994) 1747. 

[2] F.Abe, et al, CDF Collaboration, FERMILAB-PUB-94/097-E, 1994. 

[3] Z. Kunszt and F. Zwirner, Nucl. Phys. B385 (1992) 3. 

[4] B.A.Kniehl, preprint DESY 93-069, 1993; Phys. Rep. (in print). 

[5] L.R.Surguladze, preprint OITS-543 (University of Oregon, 1994); [hep-ph / 9405325| . 

[6] J.F.Gunion, H.E. Haber, G.Kane and S.Dawson, The Higgs Hunter's Guide (Addison- 
Wesley, Redwood City) 1990. 

[7] M.Drees and K.Hikasa, Phys. Rev. D41 (1990) 1547. 

[8] G.t'Hooft and M.Vehman, Nucl.Phys. B44 (1972) 189. 

[9] G.t'Hooft, Nucl.Phys. B61 (1973) 455. 
[10] W.Bardeen, A.Buras, D.Duke and T.Muta, Phys.Rev. D18 (1978) 3998. 
[11] L.R.Surguladze, FERMILAB-PUB-92/191-T, 1992. 



7 



[12] S.G.Gorishny, A.L.Kataev, S.A.Larin and L.R.Surguladze, Mod. Phys. Lett. A5 (1990) 
2703. 

[13] S.G.Gorishny, A.L.Kataev, S.A.Larin and L.R.Surguladze, Phys. Rev. D43 (1991) 1633. 
[14] D.J.Broadhurst, N.Gray and K.Schilcher, Z.Phys. C52 (1991) 111. 

[15] C.A.Dominguez and N.Pavcr, Phys.Lett. 293B (1992) 197. 

[16] J.Gasser and H.Leutwyler, Phys.Rep. 87 (1982) 77; see also S. Narison, Phys. Lett. B197 
(1987) 405. 

[17] C.A.Dominguez, C. van Gend and N.Paver, Phys.Lett. 253B (1991) 241. 
[18] J.Gasser and H.Leutwyler, Nucl.Phys. B250 (1985) 465. 
[19] T.Appelquist and J.Carazzone, Phys. Rev. Dll (1975) 2856. 



8 



